Vibrations of a Columnar Vortex in a Trapped Bose-Einstein Condensate 
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We derive a governing equation for a Kelvin wave supported on a vortex line in a Bose-Einstein 
condensate, in a rotating cylindrically symmetric parabolic trap. From this solution the Kelvin 
wave dispersion relation is determined. In the limit of an oblate trap and in the absence of longi- 
tudinal trapping our results are consistent with previous work. We show that the derived Kelvin 
wave dispersion in the general case is in quantitative agreement with numerical calculations of the 
Bogoliubov spectrum and offer a significant improvement upon previous analytical work. 
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The behavior of turbulent flows, tornadoes, mix- 
ing processes, synoptic scale weather phenomena and 
sunspots all critically depend on our understanding of 
vortex dynamics. Quantitative investigation of vortices 
started in the mid 1800's with the development of the 
Navier-Stokes equation, with the properties of a vor- 
tex being described by streamlines of vorticity. In 1880, 
Thomson (Lord Kelvin) determined the dispersion rela- 
tion [T] for a specific excitation on a vortex line. This ex- 
citation takes the form of a chiral normal mode in which 
the perturbation propagates along a vortex line and ro- 
tates about its unperturbed position, distorting the vor- 
tex line into a helical shape. 

In superfluids these vortices differ from their classical 
counterparts by having quantized circulation and a sin- 
gle line of vorticity (a vortex line) associated with them 
[2J. A quantised vortex line, in analogy with a classical 
vortex, supports Kelvin waves [2"HTT]. The first exper- 
imental investigations of Kelvin waves, in a supcrfluid, 
were in carried out in cryogenically cooled helium [2J- 
HIIH]. More recently Bose-Einstein condensates (BECs) 
have provided a new platform to investigate the proper- 
ties of quantized vortices [T2TH8] . The highly controllable 
nature of BEC systems has enabled the experimental in- 
vestigation of Kelvin waves on a single vortex line [3]. 
The behaviour of such waves has been investigated [H)l - 
[2Tj and plays a crucial role in understanding the details 
of superfluid turbulence [22H24] , In trapped systems the 
Kelvin wave dispersion for a single vortex line has been 
obtained numerically, via solving the Bogoliubov spec- 
trum for a single vortex line [2UJ 122 US]. The Kelvin 



wave dispersion relation in the limit of long wavelengths 
and in the absence of trapping is [5] 

W= ^ ln (^)' (1) 

where w is the excitation frequency of the mode, k is the 
wavenumber, M is the particle mass and r c is the vortex 
core parameter. 

A general formulation for a quantised vortex line in a 
trapped rotating BEC has proven difficult. Most meth- 
ods have relied on matched asymptotic expansions [26l - 
[29] . Koens and Martin used such a procedure to ob- 
tain a set of equations that describe the behaviour of a 
perturbed vortex line [25]. In this analysis the general 
equations [Eqs. (68) and (69) in Ref. [25J] contained an 
undetermined constant. Here we eliminate the constant 
to obtain a single equation that describes the radial posi- 
tion of a vortex line, supporting a Kelvin wave. We then 
determine the general solutions of this equation, enabling 
us to make favorable comparisons with previous results 
in limiting regimes 28 ,30 33] . From the general solution 
the Kelvin wave dispersion relation, for a single vortex 
line in a BEC in a cylindrically symmetric parabolic trap, 
is determined. This dispersion quantitatively agrees with 
numerically calculated Bogoliubov spectra [2T], in con- 
trast to previous analytic results [29"] . 

As shown in Ref. 28, 34] the equations governing the 
positional dependence of the vortex line can be obtained. 
In cylindrical coordinates, defined by the radial, p, angu- 
lar, </>, and axial, z, unit vectors, these equations are: 
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where p and describe the radial and angular coordi- 
nates of the vortex line at a given z (the dependence on 
z is implicit), Vt r is the trapping potential, I^tfI 2 is the 
Thomas-Fermi condensate density, g is the inter-particle 
interaction strength, R± is the radial Thomas-Fermi ra- 
dius, fi is the rotation frequency of the trap, and E is an 
unknown constant vector. In general, these equations re- 
late the motion of the vortex line [left hand sides (LHS) 
of Eqs. ([2| and |3])] to its distortion and coupling to 
the trapping potential [right hand sides (RHS) of Eqs. 
([2]) and Q3p] . The constant provides coupling between 
the motion of the vortex line and the trapping potential 



in the z-direction. For the equations to be dimension- 
ally consistent, E needs to take a form similar to the 
first terms on the RHS of Eqs. ^ and Hence, each 
equation has the same constant, which can be eliminated. 

In this work, we are interested in the properties of 
helical waves described by cj> — kz — Lit, in the pres- 
ence of a harmonic trapping potential Vf r (rj_ = p,z) — 
Mujj_p 2 /2 + Moj 2 z 2 /2, where lu± and u> z are the trapping 
frequencies in the r± and z directions, respectively. Ap- 
plying these conditions with p being time independent, 
the governing equation is 
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where 6 = - ln(r cV /fc 2 - d 2 p/p)/2. The LHS of Eq. @ 
defines the properties of a vortex line when p is constant. 
The first of the terms on the RHS contains the influence 
from z confinement and the other two represent the influ- 
ence the curving of the vortex line in the radial direction 
has on its motion in (f> and p respectively. 

Equation Q does not admit an analytic solution in 
a closed form. However, by assuming that: (i) d 2 p/p 
is small compared to k 2 and may be approximated by 
a constant f3 and (ii) taking the effect of ^I^tfI 2 to be 
a constant, to leading order, defined by the condensate 
chemical potential, p — g\^?TF{0, 0)| 2 , Eq. Q reduces to 
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where C = - ln(r c ^fc 2 - /?) /2 and e = 1 + 2(. In Eq. Q 
z' = z/R z , with R z being the axial Thomas-Fermi ra- 
dius, p' = p/r c , 
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where A 



m^ui' 2 ' w ^ n a being the s-wave scattering 



length, N the number of particles in the condensate and 
A = u z /w±. 



Making a change of variables such that p'(z') has the 
form f(z')^ 7 ^ Eq. ^ reduces to 
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where A = 7(5 — C)/C> H p (z) is the Hermite polyno- 
mial of order p, iFi(l;m;x) is the Kummcr confluent 
hyper-geometric function, and C\ and C2 are integration 
constants. 

Due to the axi-symmetry of the trapping potential, the 
vortex line structure described is physical when p' is not 
multivalued within the condensate radius (z = ±R Z ), 
and is symmetric or anti-symmetric across z — 0. Equa- 
tion (|8| can satisfy these conditions if C\ = 0, reducing 
the solution to the symmetric hyper-geometric function 
or if A is a positive even integer, making Ha(x) symmet- 
ric. In general there is no physical reason to restrict A 
to be an integer, hence we only consider the case where 
C\ = 0. This solution gives p' a "U" shape, centered 
around z' = 0, for example see Fig. 1(a). In this ar- 
rangement, an anti-symmetric vortex line has a higher 
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energy than its symmetric counterpart, hence this solu- 
tion does not support anti-symmetric structures. 

From Eq. ^ it is possible to check the assumption 
\(3\ < k 2 . At z — |/?| is typically of order 10~ 9 /im~ 2 
for k 2 « 10~ 2 /im~ 2 . Nevertheless /3 does depend on the 
ratio of the trapping frequencies, A, with (3 — > as A — > oo 
and j3 — > —oo as A — > 0. The divergence in (3 as A — > is 
slow, indicating that /? plays a relatively insignificant role 
in Eq. Q unless the trap is extremely prolate. In general, 
we calculate /3 self-consistently, through the definition 
P = d 2 p/p. 

Oblate Limit - Equation {[sj) simplifies in the limits of 
extremely oblate (A — > oo) and prolate (A — > 0) trapping 
potentials. In the oblate limit, 7 — > 0, £ — >• — In (r c |/c|) /2 
and £ — > 1 — In (r c |fc|), Eq. ([5| has the solution 
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where Ci and C2 are again integration constants and 
erf(a;) is the error function. Applying the symmetry and 

divergence conditions p' — (2v / C C , 2) _T + f ■ If p' is finite, 
from Eq. ^ the Kelvin wave has the frequency 
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This is consistent with previous work 30, 31 which shows 
an off centred vortex in an oblate condensate is straight 
and has a precession frequency SI — 3fiijj_ \n(R±/r c )/4:p, 
in the limit k — > 0. 

Prolate Limit - In the prolate limit (A — > 0), 7, £ and 
e all tend to 00. Hence Eq. |5]) becomes 
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Again p' only has a symmetric solution (C2 = 0). This 
causes p' to have a 1/ cos(,z) like structure similar to that 
predicted in simulations [3H [33] ■ Typical solutions of 
Eq. (12) are plotted in Fig. [l](a), demonstrating the "U" 



shape of the radial coordinate. Furthermore as — 
depends on k, this "U" shape will straighten as the mo- 
mentum of the wave increases, as seen in Fig. [lja). 
To obtain the dispersion relation for the prolate case 



Eq. ( 12 ) needs to be single- valued within the condensate. 



This gives the condition J AjC, = tt/2 or 
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Figure 1: (Colour online) (a) The vortex line shape in a pro- 
late BEC, as given by Eq. p2}, for k = 10" 2 ^m _1 (red 
dashed dotted curve), k = 10~ /im _1 (green dashed curve) 
and k = 1 /im -1 (black solid curve). We have assumed 
R ± /r c = 5, r c = 0.1 pm, ^/AJC = tt/2, CI = 0.5, and 
C2 = 0. The thick black solid line defines the Thomas-Fermi 
surface of the BEC. (b,c) The dispersion relations for Kelvin 
waves, for N = 10 4 (lower green realisations), N = 5 x 10 
(middle red realisations), and N — 5 x 10 5 (upper back reali- 
sations) Rubidium atoms in a harmonic trap: cu± — 2it x 98.5 
rad/s and lo z — 2tt x 11.8 rad/s. In (b,c) the numerical cal- 
culations (dots) [2T] and analytic results from Eq. ( |15[ ) (solid 
curves) are plotted. In (c) the dashed curves show previous 
analytic results For the solid curves and dashed curves 

in (b,c) the following scaling parameters have been used for 
r c and k: a = 0.134, 0.257, and 0.264 for N = 5 x 10 5 , 5 x 10 4 , 
and 10 4 , respectively, and s — (XRz/l^^ 1 ^ 2 for all N. 



which is an approximation to the dispersion relation of 
helical waves in a prolate trap, and is distinct from the 
solution for uj z = 0. This occurs because any trapping in 
z provides a scale over which the vortex line can curve. 
Hence, when there is no trapping in z, e/7 = (/"f = 0, 
turning Eq. ^ into 7p' 2 = 0. Therefore, for the case of 
lj z = the dispersion relation is the same as the oblate 



limit, given by Eq. (10), i.e. the vortex line is straight 

General Solution - To determine the general disper- 
sion relation from Eq. (JsJ) a similar treatment to the ex- 
tremely prolate trap is employed. As the p' solution is 
a hyper-geometric function to a negative power, the ze- 
ros of the hyper-geometric function must lie outside the 
condensate. The zeros of the confluent hyper-geometric 
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function iFi(l;m;x) can be approximated by 
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where Xq is the approximate x value of the rth zero. 



Combining Eq. ( 14 ) with the form of the hyper-geometric 
function in Eq. and stipulating that the first zero is 
at z' = 1 the dispersion relation becomes 



W2D 



8 M (1 + C) 



(15) 



or equivalently, for (3/k 2 —> 0, 



UJq + U3\ + 



fr fc21n (A) 

2M 



n, 



(16) 



where 



3huj 



w 



4a*, 

fiw 2 ln(r c |fc|) [4 + 7T 2 ln(r c |fc|)] 
16/i [2 -In (r c |fc|)] 



Equations (15 1 and (16) indicate that helical waves, 



in a parabolic trap, obey the usual dispersion relation 
[Eq. Q], with two constants wo, from confinement in p, 
and wi, from confinement in z. The form of luq matches 



that for extremely oblate traps [Eq. (10)], while ui\ con- 
tains new behaviour. Essentially, u>i is constant, with 
weak logarithmic dependence on k, and becomes larger 
increases. As A — > oo Eq. 



as oj z 



(15 



does not replicate 
the dispersion in the oblate trapping limit. This is be- 



cause in the very oblate limit Eq. ( 14 ) fails to predict the 
location of the zeros accurately. 

To test the validity of the above analysis we now com- 



pare the solutions of Eq. ( 15 ) with numerical calcula- 
tions [21 J. The dispersion relation of Kelvin waves in 
a prolate condensate for different particle numbers cor- 
responding to numerical calculations |21| are shown in 
Fig. [ljb,c) [dots]. Previous analytic predictions in this 
limit, see Eq. (70) in Ref. [29], poorly replicated theses 
results [dashed curves in Fig. Wc)]. To compare these 
numerical results with Eq. |l5]), r c needs to first be con- 
sidered. 

The core parameter, r c , implicit in Eq. (15), char- 
acterizes the vortex core size. In a trapped BEC the 
healing length, £, is of order the vortex core radius, 
which is a function of position. As such we define the 
core parameter r c = at; to be some fraction, a, of the 
healing length averaged over the Thomas-Fermi volume: 
£ = (R\R z /QaN) 1 / 2 . To compare with numerical results 
we allow a to be a free parameter, of order 1. 

In Figs. l(b,c) (solid curves) we plot the Kelvin wave 
dispersion, Eq. (15), where we have rescaled k — > sk — 



k s . Due the inhomogeneity of the condensate we expect 
k to vary spatially. As such this rescaling is motivated 
by the observation [21 that in numerical calculations k 
varies along the vortex line. Interestingly, we find that 
the matching between analytical results and numerical 
calculations is optimized for s = (A-Rz/Z^) 1 / 2 , where 
l Uz = \J h/(muj z ) is the harmonic oscillator length in z. 

Figure. [l](c) shows that Eq. ( 15 1 (solid curves) is a vast 
improvement compared with the previous analytic pre- 
dictions (dashed curves). Additionally, Fig. [ljb) shows 
excellent quantitative agreement between the numerical 
and analytical results, with the agreement improving as 
the number of particles increases and the BEC becomes 
more Thomas-Fermi like. The general features of the 
dispersion relation are dominated by: (i) a frequency 
shift and (ii) a functional form oc —k 2 \nr c k. The shift 
is dominated by the confinement in p, given by wo, with 
a small contribution from the confinement in z, given 
by the second term in Eq. (15). The functional depen- 



dence of lj is essentially dominated by the solution for 
a vortex line in an un-trapped condensate, Eq. ([lj, also 
with a small correction arising from the confinement in 



z, given by the second term in Eq. (15). We note that 



the self-consistent determination of j3 only influences the 
dispersion for k < 10~ 4 praT 1 and hence, in Figs. [ljb,c), 
plays no role over the wave-vectors considered. 

In summary, we have analytically determined the 
Kelvin wave dispersion relation for a vortex line in a 
BEC trapped in a cylindrically symmetric parabolic trap, 
Eq. (15). This result quantitatively agrees with numer- 
ical calculations, in stark contrast to previous descrip- 
tions. We also find that the dispersion relation in the 
oblate trapping limit, Eq. (10), coincides with previous 



analytical work. In general these results are derived from 
the governing equation of motion for Kelvin waves on a 
quantized vortex line, Eq. ([5]). Equation ^ successfully 
predicts the behaviour in very oblate and prolate con- 
densates, and mathematically shows the link between the 
vortex line shape with no trapping in z and strong trap- 
ping in z. Within the context of experimental activity 
in the study of excitations of quantized vortex lines in 
trapped BECs HH |36l |37j , this work provides a sim- 
ple analytic tool for the analysis of Kelvin waves. Ad- 
ditionally, given the close agreement between numerical 



calculations and Eq. (15) we expect future experimental 



measurements of the Kelvin wave spectra to be in close 
agreement with our analytical description. 
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